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ABSTRACT 

I derive analytic scalings for coronagraphic imaging searches for extrasolar planets. 
I compute the efficiency of detecting planets about any given star, and from this 
compute dimensionless distribution functions for the detected planets as a function of 
planet-star distance and distance to the host stars. I find the following for blind planet 
surveys: (1) the optimum wavelength is between 4000-5000 A for Earth-like planets 
and 4200-5800 A for Jovian planets; (2) between 21-32% of the number of planets per 
decade of radius can be detected with an optimized survey; (3) target stars should be 
ranked from greatest to least by their luminosity divided by distance to the sixth or 
eighth power, depending on the dominant source of noise for the survey; (4) surveys 
targeting all main sequence stars will detect ~3 times as many planets as surveys only 
targeting G-type stars; and (5) stellar populations with different metallicities should 
have exposure times that vary with the cube of the metallicity. I apply these results 
to the current suite of proposed coronagraphic satellite telescopes, of which TPF-C 
is the most powerful, but a much smaller telescope, TOPS, may have a significant 
chance of detecting Earth-sized planets due to its small inner working angle and high 
throughput. The most significant uncertainty in these results is the noise contribution 
of Exo-zodiacal light. These results can be applied to designing coronagraphs, com- 
paring proposed telescope designs, optimizing the observing strategies, determining 
the properties of detected planet populations, and selecting target stars. 
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1 INTRODUCTION 

Now that extrasolar p l anets have been discovered with various indirect means (radial velocity, transits and microlensing, 
IMavor fc Quelozlll99.4 IKonacki et alJ 1200.4 IBond fc et alJ l2004h . it is now a major goal of astronomy to directly detect 
extrasolar planets. Two mea ns of imaging planets have been proposed: optical c o ronagraphy to cap t ure light from the 
star scattered by the planet l|Nisenson fc Papalioliosl l200lt iKuchner fc Traublbooi Easdin et aJ 1200.4 iBrown et all 1200.4 
Stapelfeldt et alJl200Sh and infrared interferometry to detect thermal emission from the planet 1 jLeeer et alJll996l : iFridlundl 
2004; Kal teneeeer fc Frid lund 2005J). Here I carry out a detailed analytic analysis of the prospects for planet detection using 
optical coronagraphy with space-based telescopes. 

A coronagraph enables planet detection by suppressing the light from the central star and suppressing the wings of the 
point spread function. The remaining stellar light is due to diffraction and scattering due to imperfections in the optics, which 
in some designs can be suppressed further by differencing images taken at different roll angles so that the speckles can be 
subtracted while the planet remains. Typical coronagraphic designs have an "inner working angle" within which the stellar 
PSF is too bright and the planet light is suppressed by the coronagraph. This sets the minimum angular scale at which 
planets can be detected relative to their host stars. Secondly, a coronagraph has a contrast between the wings of the PSF 
and the unobscured central brightness of the star, which, in addition to flux from local zodiacal light in our solar system and 
"exo-zodiacal" emission due to scattering by dust in the observed planetary system, sets the background at the position of 



rhe recent detections of the eclipse of planet b y their host stars also constitute direct detect ic 
from the star and which are from the planet (Charbonncau et al. 2005; Dcmine et al. 2005). 
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the planet. These backgrounds set the maximum planet-star separation at which the planet can be detected as the planet 
flux decreases as the separation increases. Since the inner working angle is proportional to the wavelength of light, at shorter 
wavelengths planets can be detected closer to their host stars (where they are brighter) or around host stars at greater 
distances from Earth giving more stars that may be probed. However, the albedo and the stellar flux eventually decrease 
at short enough wavelengths, so at some point it does not pay to have higher angular resolution. The PS F contrast rises 
as th e inverse of the square of the wavelength (due to the increase in phase error for shorter wavelengths, lAneel fc Wool] 
1998), which partly counteracts the benefits of short wavelengths if the wings of the PSF dominate the noise. If zodiacal or 
exo-zodiacal light dominate the noise, then shorter wavelengths means that the planet PSF will be more concentrated and 
thus less affected by the zodi/exo-zodi surface brightness. To compute the optimum wavelength of observation requires a 
knowledge of the sample stars, coronagraph properties, zodiacal and exo-zodicacal light, and survey parameters, as well as a 
model of the noise and coronagraph response. It is the goal of this paper to derive these quantities in as simple a manner as 
possible to elucidate the optimization of coronagraphic surveys for extrasolar planets. 

As a measure of the success of a coronagraphic survey I use a straightforward metric: the total number of planets detected 
with specified properties. This allows a quantitative comparison of various telescope designs, survey strategies, observation 
wavelength, and choice of stellar targets. For the purposes of this study, I confine the analysis to blind, single-visit surveys, 
i.e. for which there are no known planets and each star is observed a single time. 

Computations similar to those carried out here have been presented bv lBrownl I^OOfiFl . the primary difference being that 
Brown has focused on habitable terrestrial planets and his computations are numerical (using Monte Carlo realizations of 
simulated surveys) which require great computational expense to explore parameter space for different surveys. The formulae 
presented here target detection of planets at any separation from their host star and are analytic, allowing a rapid exploration 
of parameter space and faster optimization of surveys. 

After introducing my assumptions about the planets (fj^J and coronagraphs (Q, I compute the number of stars that can 
be surveyed (©, the efficiency of detecting planets about a given star O, and the planet and star distance distributions 
(101. I show how to choose the (dimensionless) survey volume ([JHJ, how to rank target stars (©, an d how to optimize the 
wavelength of observation f m.00 . I use my results to compare various proposed coronagraphic imaging satellites ( jfljl . I end 
with a discussion of how to weight observing time by metallicity ( 812(1 . In 8131 1 discuss the results and outline topics for 
further study. Two tables in the appendix summarize the notation used in the paper. 



2 PLANET DETECTION REQUIREMENTS 

To detect a planet with a coronagraphic telescope requires 

• sufficient signal-to-noise, S/N ^ (S/N)det- This should be chosen to be large enough that one is confident that a source 
is not due to a statistical fluctuation, so (S/N) de t ~ 5 — 10 should be sufficient given that tens to th ousands of stars wi ll be 
surveyed (although higher signal-to-noise may be required to provide convincing evidence of a planet. iGould et alJ l2004V 

• sufficient angular resolution, 9 — rsina/D ^ 8iwa, where D is the distance to the star, 9 is the angular separation on 
the sky of the planet and star, r is the planet-star separation, ^ a ^ n is the phase angle of the planet (a = at full phase), 
and 9iwa is the inner-working-angle of the telescope within which a planet cannot be detected near the star. 

These two conditions lead to the differential number of planets a survey can detect, 

''- V " =\^sD*^%-g-H{9-9 1W A)H{S/N-S/N det ), (1) 



drdR p dadMdXdD 2 dM dX drdR p 

where Q, s is the survey solid-angle, D is the distance to surveyed stars (Q S D 2 is the volume surveyed), dn/dM is the stellar 
mass function (I assume in this paper that host stars are on the main-sequence), dp/dX is the stellar metallicity probability 
distribution (Jf=[Fe/H]), df / (drdR p ) (defined below) is the frequency distribution of planets as a function of planet-star 
separation r and planet radius R p , and H{x) is the Heaviside step function. I have assumed azimuthal symmetry for the 
detectability of planets around a given star, so the angular detection limit just depends on 27rr 2 sin a/(47rr 2 ). I have assumed 
that the stellar metallicity distribution is independent of stellar mass, hence the separate dependence on 4?j and -jfe, and I 
have assumed that the planet size distribution is independent of stellar properties. The two step functions enforce 6 > 9iwa 
and S/N > S/Ndet- In the next section I specify my assumptions about the planet properties used in integrating this equation. 



3 PLANET PARAMETER PRESUMPTIONS 

In order to make a straightforward study of the dependence of coronagraphic searches on the properties of the coronagraph 
and telescope, I make some simple assumptions about the planet distribution and physical properties of the reflecting planets. 



© 0000 RAS, MNRAS 000, 000-000 



Planet Searches 3 



I assume that the distribution of planets is independent of stellar spectral type, that the survey targets planets of a particular 
size R p , and that the planet frequency is constant with the logarithm of radius, 

df _ ho 



dlnrdR'p In 10 



s(K - R P ), (2) 



where /io is average number of planets per decade of radius and 5(x) is the Dirac delta function. Since I am assuming single- 
visit observations of each star, I can ignore the orbital parameters of each planet (such as eccentricity) and simply treat the 
density of planets near each star as a power law in radius once averaged over planet orbital elements. This implies a spatial 
density of n(r) = /io/(47rr 3 In 10). Now, if the planets are distributed between an inner radius, n n , and an outer radius, r ou t, 
then the total fraction of stars with planets is f p = /io[ln(r out /ri„)/ ln(10)]. I make this choice of planet distribution since the 
distribution of giant extrasolar planet periods is nearly proportional to P- 1 iStepinski fc BlacklEoOlt iTabachnik fc Tremainei 
l200llKuchneill2004) whi ch implies a constant number of planets per logarithmic radius interval, although one study suggests 
a shallower dependence jLineweaver fc Gretheill2003h . We currently have no information on the distribution of terrestrial 
planets, but the constant In r distribution seems to be a sensible choice as it is not biased towards detecting planets at either 
larger or smaller radius. 

I compute the planet brightness with several assumptions. First, I assume that the brightness of the planet scales as the 
inverse square of its distance from the star. Second, I assume that the geometric albedo has a fixed spectral shape independent 
of planet-star separation. Third, I assume that the phase function is independent of planet-star separation and wavelength. I 
assume that all of these properties are independent of the planet mass, stellar spectral type, age and metallicity. Then, for a 
single exposure of a star with a planet, the number of photons detected from the planet is 

Q p ~Q* Px (?A\(a), (3) 

where p\ is the geometric albedo (the fraction of the star's flux reflected at full phase) of the planet at the observed wavelength 
A , $(a) is the phase function (defined to be 1 at full phase when the star-planet-observer angle a = 0), and Q* is the number 
of photons from the star. 

Planetary phase functions are fairly complex, but given the generality of this paper, I use the phase function <!>(«) = 
cos 4 (a/2). I refer to this as the "quasi-Lambert" phase function as it approximates the Lambert phase function, $(a) = (sina+ 
(n — a) cosa)/-7r, which is the phase function of diffuse scattering. The quasi-Lambert phase function has a mathematically 
convenient form that allows for analytic solution of the planet detection efficiency. 



4 CORONAGRAPH AND SURVEY ASSUMPTIONS 

I assume that the coronagraphic search is carried out with direct imaging at a central wavelength A with a fractional band- 
pass A In A < 1 with a mean throughput of e(A). The parameter e includes all inefficiencies due to, e.g., partial obscuration 
of the mirror by the secondary, partial reflectance, partial coronagraphic obscuration, quantum inefficiency, observational 
overheads due to target acquisition and slewing, readout time, and multiple exposures for covering multiple angular regions 
for non-axisymmetric PSFs. I assume a circular telescope aperture of radius Rtei with an inner working angle 

Oiwa = e>iwA\/(2Rtei), (4) 

where Oiwa is a dimensionless number, typically 1.5-5, which relates the inner working angle and the ratio of the wavelength 
to the telescope diameter. 

I assume that the telescope can achieve a contrast, C(A), which is constant in an annular region Oiwa < < Bow a, 
where Oowa 3> Oiwa is the outer working angle. Th e contrast is define d as the intensity ratio at a point in the PSF to the 
intensity of the centre of an unocculted stellar PSF jKasdin et aj feoOS). Since the ph ase errors scale as A , I assume that 
the contrast ratio scales as C(A) = Cq(Aq /A) 2 , where Ao is som e reference wavelength JlVlalbet et alJflQQfl^ . This dependence 
results from gaussian errors in the phase ijAneel fc WooHlll998^ . 

I assume that the telescope carries out a blind survey of nearby stars (i.e. not choosing stars that are known to have 
planets). I assume a constant exposure time of each star, T exp (if each planet is observed multiple times at different roll angles 
for speckle subtraction or covering the detection zone of the PSF, then T exp is the total time for exposures which include the 
planet). I assume that the survey has a fixed duration, T s , which is typically of order several years, so that the total number 
of stars surveyed is 

N s = T s /T exp (5) 

(neglecting overheads, which are likely small due to long exposure times). Finally, I assume that each star is observed at one 
epoch so that the planet is at a fixed position. I relax the constant exposure time assumption later. 



© 0000 RAS, MNRAS 000, 000-000 



4 Agol 



5 SIGNAL-TO-NOISE RATIO 

The total signal-to-noise for an observation scales as 



S/N = T7 —„ ^ — ^ (6) 



Qv 

(Qz + Qez + Qpsf + Qb) 1 

where Q p is the total number of detected photons from the planet, Qz is the noise contribution from zodiacal light, Qez is 
the noise contribution from exo-zodiacal light, Qpsf is the noise contribution from the wings of the stellar PSF, and Qb is 
the noise contribution from other sources (e.g. read-noise and dark current noise). This equation assumes the large-number 
photon limit applies, and since I am computing the number of planets that can be detected above a given signal-to-noise, I 
neglect the noise due to the planet (for planet characterization this should be included). For the remainder of this paper I 
will ignore Qb as I expect that detectors will be designed so that instrumental noise does not dominate. 
The total number of photons detected from the unobscured star is 

Q* 7ri& ! e i |^A In \h~ 1 T exp (7) 

where L v is the star's specific luminosity (in erg s~ Hz -1 ), h is Planck's constant and the other quantities are defined above. 
The noise due to zodiacal light is given by 

Qz = ^(lAU) 2 r z b fac (X/D tel ) , (8) 

where tz (units of sr _1 ) relates the solar flux to the zodiacal light surface brightness in the direction of ecliptic longitude and 
latitude (A ec ,/3ec) and Sf ac measures the PSF sharpness in units of \/D te i (e.g. for optimal PSF fitting with an Airy disc 
point-spread-function, Sf ac — 2.1 ). The noise due to exo-zodiacal light is given by 

,o irRt el eL v Ala\h 1 T exp /\/r> \2 / n \ 

Qez = 7—7 : 75 TEZOf ac {\ /Utel) , (9) 

47r(rsma) 2 

where tez (units of sr ) relates the surface brightness of the exo-zodiacal light at the position of the planet to the flux of 
the star at the tangent point (integrated along the line of sight). Finally, the noise due to the wings of the PSF is given by 

Qpsf = Q*C(\)S fac - (10) 



6 SIMPLE SURVEY SCALINGS 

Having outlined my various survey assumptions, I now derive some scaling parameters for coronagraphic surveys. I first assume 
that the noise is dominated either by Qz, Qez, or Qpsf', in >I6.5I I will combine these to determine the scalings when all three 
sources of noise contribute. 

There is a toroidal-shaped region surrounding the star within which the planet satisfies the signal-to-noise detection 
criterion (Figure^. As the star becomes more distant the signal-to-noise for planets at separation r decreases (oc I)! -1 ' -2 ' -1 / 2 ' 
in the [Zodi, Exo-zodi, PSF] limited cases), so the signal-to- noise detection region shrinks for a fixed exposure time. The angular 
resolution limit corresponds to a cylinder surrounding the star with the axis along the line of sight with a radius which grows 
as distance set by the inner working angle. Inside the torus and outside the cylinder is the planet detection zone. At some 
distance, D max (T exp ) these two limits meet so that the torus is contained entirely within the cylinder - no planets can be 
detected for stars at D > D ma i(T e ip) for a particular set of planet, stellar and telescope parameters. The distance D ma x sets 
the scale for the volume of stars that can be surveyed of a particular stellar type, so I next compute this for each of the three 
noise limits. 



6.1 Zodiacal dominated noise 

When Qz S> Qez, Qpsf, to detect a planet requires Q p /Qz ^ (S/N)d e t- This implies 



r < 



The resolution limit requires 



1/4 



= r 2 ,z- 



i (sin a) 



A 
~Dtei 



D 



ri. 



(11) 



(12) 



Since detection requires n ^ r ^ r^.z, the limiting distance D max is set by r\ = r 2 ,z at a max = 60° which is where the planet 
has the largest angular separation from the star (this is where the flux-limit and resolution limit meet in panel d of Figure 
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Figure 1. Cylindrical cross section of the region where planets can be detected for stars at a distance D from the Sun for (a) D = 
0.25D max ; (b) D = 0.5D max ; (b) D = 0.75D max ; (d) D = D max . Each panel is a 2-D cut through a 3-D region surrounding the 
star at the origin of each panel, where the units of the axes are Qiw A^max- The observer is at [x, z) = (0, —DOJ^ A D^a X ), and the 
inner working distance is indicated with vertical solid lines. The solid lines are for the Zodical dust noise limit, the dashed lines for the 
Exo-zodiacal noise limit, and the dotted lines for the PSF-limited noise limit. Note that for each noise limit, D max will be different, so 
each panel compares the three noise limits at different physical distances. The shaded regions are where a planet is detectable outside 
the inner working angle and inside the flux limit (I have only shaded one of the noise limits in each panel for clarity). Since panel (d) is 



at D = D„ 



no planets can be detected. 



0. This gives 



where 5 = 1/8 and 



D max — Dz = Dte 



LlT exp (lAU) 2 



\ 6 L e ,„TZ@ 4 n 



Q = -^Rppl$(a m a X ) 2 sin 4 a max eAln\h (S/N) d 2 t S fac 



(13) 
(14) 



6.2 Exo-zodiacal dominated noise 

The noise due to exo-zodiacal light is difficult to estimate as the zodiacal dust content of other old stellar systems is un- 
known, except for a ha ndful of infrared constraints with a sensitivity to dust ~75 times as bright as our zodiacal dust cloud 
l|Beichman et al.ll200tj) . So, the largest uncertainty in planet detection efficiency is due to exo-zodiacal dust; until better 
constraints become available we will have to rely on educated guesses to estimate the contribution of this to planet detection 
noise. 

For a planetary system with a dust cloud similar to that of the solar system, the surface brightness declines steeply with 
distance from the st ar, oc r~ 23 , and var ies significantly with the inclination angle and azimuth of the disk due to flatness of 
zodiacal dust cloud llLeinert et aljfl998i) . I will consider the case of largest noise: an edge-on disk. In this limit the equations 
simplify so that the exo-zodiacal surface brightness is just a function of the angular separation of the planet from the star. 
The detection limit for signal-to-noise, Qp/QgJ ^ (S/N)d e t, translates into a similar limit on radius which is 

r < [R 4 p pl$(a) 2 Q*D- 2 sin 2 a T E lSjlD 2 tel \- 2 (S/N)^ 2 t ] 1/2 = r 2 , EZ . (15) 
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This equation applies when tez is independent of r (which is nearly true for our zodiacal cloud). The resolution limit is the 
same, so setting r\ = T2.ez at a max (by coincidence, a max = 60° in this case as well) gives the maximum distance 

1 s 



(16) 



D max = Dpsf = D t . 



(18) 



D max — DeZ = Dtel 

' ~'~'IWA. 

where 8 — 1/6 here. For some stars with zodiacal dust that is more face-on, Dez will be greater than this estimate. However, 
given the uncertainties in estimating tez, a more complete calculation is unwarranted at this point. 

6.3 PSF dominated noise 

The signal-to-noise limit in the PSF noise dominated case, Qp/Qp^p ^ (S/N)det, leads to the requirement that 

r < [Q^C-'SJ^plR^ia^S/N)^} 1/4 = r 2 , PS F (17) 
Setting n = T2.PSF at a max allows us to solve for 

T T 1 

where 5 = 1/6 here. With these expressions for D max in three limits, I turn the planet detection expression into a dimensionless 
equation in f3 

6.4 Maximum number of stars 

I now compute the the maximum number of stars that may be surveyed of a particular spectral type. If all stars surveyed are 
of a specific spectral type, i.e. dn/dM' = n*5(M' — M), then the total number of stars surveyed is N s = Q s n„D^ /3, where 
n* is the number density of stars of the given spectral type, assumed to be constant, and D s < D max is the limiting distance 
of surveyed stars. The total number of stars surveyed is also a function of the survey duration, N s = T s /T exp . 

The maximum number of stars that can be surveyed, N max , can be found by equating these two expressions for N„ and 
setting D s — D max — D[ Z ,ez,psf\- I can solve for the maximum number of stars that can be surveyed in the three different 
noise limits 

N max = N [Zt Ez,psF]{T s ) = [\£lsn*D {ZiEZ)PSF] {T s f} , (19) 

where D[ Z ,ez,psf]{Ts) means that T exp is replaced by T s in equations l|130 . I|16|l . or l|18|l . 

Since N max sets the characteristic scale for the number of stars to survey for a particular set of telescope, stellar, and 
planet parameters as well as total survey duration, it makes sense to scale the actual number of stars surveyed as 

Af = N 3 /N max , (20) 

where N < 1. In 01 determine the value of M which maximizes the number of detected planets and I estimate the total 
number of planets that can be detected, Ndet, but it should be clear that Ndet scales with the maximum number of stars that 
can be surveyed, Ndet oc N max , since the more stars that can be surveyed, the more planets can be detected. Although N max 
is defined only for a specific spectral type, in JHll generalize to surveys of stars with multiple spectral types. 



6.5 Summary 

When all three sources of noise contribute, the general equation for D max (T exp ) becomes: 

D 8 D 6 D 6 

^ m ax | ^ m ax , ^ ma x 



I max , ^ max -i ( r )~\\ 



n8 n6 1 n6 

^z ^EZ ^PSF 



which is a quartic equation in D max with solution 



1/2 



S = ^ (VV 4 + 256/27 -V- 

D% (D% z +D%sf) (99 , 

1J EZ 1 -'PSF 



The ratio of the maximum distances in the three noise limits are 



T> ^ Z _ ylrnl/4p.-l/6 x -l/12 1/6^-1/8 r -l/24 r ,-l/24 T -l/24 f -1/8 
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v 2 = = (iAc/) 1 / 4 e^ A A- s / 1 V^v^y 12 g-i/ 2 4 r -i /2 4 C7 i/ 6A V3 L -i/8. (23) 

UpSF 

The weak dependence on L„ in these equations means that different spectrai types wili typicaiiy have simiiar noise properties. 

Typicaiiy more than one source of noise dominates. However, the fiux iimits scaie as r 2 oc Dl -1 ' -2 - -1 / 2 ! in the Zodi, 
Exo-Zodi, and PSF dominated cases, respectiveiy, so if two noise limits are comparable near D max , whichever grows fastest 
will dominate at smaller D. 



* m a x 



N \ 11/3 / N \ 3 f N \ 3 

+ luTT" + at = L - \ zi > 



N z ) \ N EZ J \NpsfJ 

1 /3 

This equation requires numerical solution of an eleventh order polynomial for NrAax ; however, to an excellent approximation, 

-3/10 



N max = (Nz 10/3 +N-™ /3 +N-ir) ■ (25) 



7 DIMENSIONLESS DISTANCE DISTRIBUTIONS 

The equations describing the detectable planets can be made dimensionless by scaling the star distance, D, as 7/ = D/D max (T exp ) 
and the planet-star separation, r, as £ = r /r max where 

_ OlWADrnax 

' max — . ■ y^jKJ J 

SHI O^max 

Integrating the planet number density, df/dr over the volume contained between the flux limits and inner working angle gives 
an expectation value for the number of planets detected about a star at a distance r\ of 



dN det _ ^ ^ dn dp 



2w ( sinada [ dr^-H{6 - 6 IWA )H{S / N - {S / N) det ) 
Jo Jo dr 



dMdXdD dM dX 

where E{if) is the expectation value of the number of planets of radius _R P detected for an observation of a single star with 
mass M, metallicity X, and distance D — rjD rnax (T exp ). 



7.1 Expectation value as a function of distance 

In the limit that one source of noise dominates at all distances, I can compute the integral of E(rj) over a and r, finding 

E{v) = ^/lo(lnlO)- 1 [3 In 

(1 + A) [l±(2A- 1 -l) 1/2 ] , 
3s 3r x 1/2 



x± = 



r = [ s (i + vT^7)] 1/3 , (28) 

for 0^7/^1, where s = r/ 1 ^ 26 ^ and = 1/2 applies in the PSF and zodi limits, while tj> = 1 in the exo-zodi limit. This 
can be approximated accurately by E{rj) = 0/io 61 ^ l 3 10 s~ 1 ^ 6 (l — s 1//3 ) 3//2 . In carrying out this integral, I have neglected the 
cutoffs in the planet distribution at r out and ri„ which is appropriate if the peak of the detected planet distribution is well 
within these radii. If df/dlnr is a power law with radius rather than a constant, then E(rf) can be expressed in terms of 
Hypergeometric functions; however, for the reasons discussed above I stick to using a uniform lnr distribution. These three 
functions are plotted in Figure [5] which both show a rise at small r\ as a larger volume is probed and fall at r\ < 1 as the 
detection threshold is reached. 

In the limit J\f = 1, the total planet fraction can be computed analytically for the three noise limits. The fraction of 
surveyed stars with detected planets is 

!57r2- 7/2 3~ 5/4 HlO)]" 1 ~ 0.153 Zodi-limit 
H [ln(lO)]" 1 ~ 0.206 Exo-zodi limit (29) 

^ [ln(10)] _1 ~ 0.103 PSF-limit 

So, depending on which noise-limit dominates 10-20% of planets per decade of radius will be detected around the stars that 
are surveyed. The detection inefficiency is due to both the spatial resolution limit and the flux limit which allow detection of 
planets within a limited range of radii, which I compute in the next section. 



N det 

N< max 
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Figure 2. The expectation value for the number of planets detected as a function of distance for Zodi-limitcd (solid), Exo-zodi limited 
(dashed), and PSF-limited (dotted) surveys multiplied by the volume factor Airr) 2 as a function of rj = D/D max . 



In the general case when all sources of noise contribute I have to integrate the dimensionless detection equations numer- 
ically. Denning 



ijz 



r)EZ = 



rjpsF 



~Dz~' 

Dmax 

Dez ' 



DpSF 

where < t]z,1]ez,tipsf < 1 since D max < D z ,D E z,D P sf, then, 

8,6 6 -, 

Vz + Vez + Vpsf = 1, 
which follows from equation 12 U . The equations for ri and r-z become 



(30) 



(31) 



r-2 



sm a 

TlEZ + ^~\Vz + VpSfV-' 2 ) ( #T 



1/2 



~1 1/2 



2 (r?I + ^ 2 4sf) 



(32) 



which is just a way of rewriting the resolution and signal-to-noise limits on radius. Using equation 1311 I can eliminate rjz 
so that r2/r m ax is just a function of rj, t]ez, rjpsF and a. With these definitions, the planet expectation value for a star at 
distance rj is given by 



E(t],t]ez,tipsf) = 



2 In 10 



r-2 



x + /2-l 

d(cos a) In I — I , 

_/2-l V r i, 



(33) 



where sc± = 4cos 2 ^ are given by equation (I28L but with s = ?7 3 [l — ??|(1 — ^ 2 )] ^ 2 (a± are the angles where the resolution 
limit crosses the flux limit). I have not found an analytic solution to this integral, so I integrate it numerically below. The 
shape is qualitatively the same as in the three noise limits since all three have very similar shapes. 
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1.000 f 




Figure 3. The expected distribution of the detected planet-star separations £ = r sin Qmax/[SiwADma.x] for the Zodi-limit (solid line), 
Exo-zodi limit (dashed line), and PSF-limit (dotted line) for surveys of a single stellar type, constant exposure time per star, and jV = 1. 



7.2 Radius distribution of planets 



In the limit Af = 1, the radial distribution of detected planets can be computed analytically in the Zodi and PSF noise limited 
cases, assuming that all stars observed are of a particular spectral type (e.g. G stars). For the Zodi limit, 



dN d 



£l s Tl*D rnax f\\ 

3 In 10 



C 4 (i-C 4 )(i + 9C 4 ) , C 2 - 

(1 + 3C 4 ) 3 2^3 



3< 4 



3C 4 + 1 



(34) 



The Exo-zodi case can be computed analytically in terms of Hypergeometric functions, but an excellent analytic approximation 
is 

dN det _ n 2 n a n«D 3 max f 10 ( 1 \ 9 
d( 4 In 10 U _2/9 + C 7/ V ' 

where g = 1.9 provides a good fit to the numerical curve. For the PSF limit, 



dN det 
dC 

w 
£ 



2( In 10 

C 3 [s - s- 1 ] , 



w" (1 3w\ vf ( 3w 



1 - 6w + 



9w* 



+ 



27 



3w 



1/3 



(36) 



These distributions are plotted in Figure El All three distributions peak near £ = 1 since the largest number of stars surveyed 
are near D ma x where the sensitivity is concentrated near f = 1 since this is near where the resolution and flux limits meet. 
Each of these curves approaches the same value at small f since the planets closest to the star are limited by resolution, which 
is independent of the phase function, scaling as oc £ 2 . At large f, planet detection is flux-limited which scales as T2 oc D~ b 
where b = [1, 2, 1/2] in the Zodi, Exo-zodi, and PSF noise limits. The planets at largest separation are detected for the closest 
stars, so it is straightforward to show that dN det /d(, oc £-( 3 / 6+1 ) = ^I- 4 .- 5 / 2 -- 7 ! m the Zodi, PSF and EZ limits, which is 
indeed the slope of the lines at large £• Since the EZ flux limit grows most slowly with radius, the observed slope will generally 
be oc C" 7 . 



8 MAXIMIZING OVER EXPOSURE TIME AND SURVEY VOLUME 

I assume for the time being that the telescope observes a single type of star at a single wavelength. Then, if fewer stars are 
observed out to a distance D s < D max (i.e. N < 1), more time can be spent on each star making it possible to detect more 
planets since these stars are closer. Furthermore, it may pay to vary the exposure time with the distance to the star, say 
T exp oc D~ K . Well, it turns out that choosing k = and N = 1 is close to optimal. By increasing « to ~0.5 (i.e. spending 
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more time on closer stars) the total number of planets found can be increased by only ~0.5% compared to the K = case 
(for either S — 1/6 or 1/8). The reason is that the most planets are detected near D s due to the larger volume at larger 
distance, but because these stars are more numerous, they also primarily determine what the exposure time should be, so 
the total planet sensitivity at that distance is roughly constant. By making k slightly negative, then you can detect slightly 
more planets around closer stars, but this is counterbalanced by the fewer planets detected around more distant stars, and 
thus the total number of planets detected remains nearly constant. One advantage of spending more time on closer stars is 
that there will be more detections at higher signal-to-noise. Specifically, the number of stars with S/N above (S/N)det scales 
as (S/N)^ S ■ Since this has a weak dependence on k, the gain is not significant, so for the rest of the paper I assume a 
uniform exposure time with distance, i.e. k = 0. 

A reduction in the number of stars can lead to an improvement in the number of planets detected as more time is spent 
observing closer stars. Figure [I] shows the dependence of Ndet on M computed from 

= 3AT 3 * / d VV 2 E(r,) (37) 



N„ 

for each of the three noise limits. For the zodi-limited case, I find M = 0.59 optimizes the total number of planets detected 
with an increase in the number of planets by ~11% over the value at N = 1, while in the Exo-zodi and PSF limited cases, 
N — 0.54 is optimum with an increase in the number of detected planets by ~16% over the value at M — 1. Even though 
only about half the number of stars are surveyed, the planet detection fraction roughly doubles causing a larger number of 
planets detected. In terms of distance, this corresponds to observing stars out to D s ~ 0.75D max rather than out to D max . 
For the general case where all three sources of noise contribute, one must optimize 

D 3 (KT\ f^ 1 " D »«..(i)/D™-M 
N det {M,Vt,V 2 ) = (n s n*Dl ax (l)) j d VV 2 E( V ,VEz,VPSF), (38) 

where I have evaluated 2?i,2?2 at A/" = 1, i.e., at T exp = T s /N max . The first term in parentheses is equal to 2>N max . The 
Zodi-noise dominated case is for T>\ = T>2 = or r\z — 1,T)ez = 0,i]psf = 0, while the Exo-zodi noise dominated case is 
for T>i — oo, T>2 = or r\z = 0,tiez ~ 1,T)psf — 0, and the PSF-noise dominated case is for T>\ — and T>2 — oo or 
rjz ~ Q,TjEZ = Q,TjpsF = 1. The quantities (r)EZ,T]psF) can be expressed as a function of [A/\ T>\, Di\ with the relations 

Vez = y 



Vpsf = y T-> 



2 ■ 



V = I 1 - Vez - Vpsf) 
y 4, +Ky 3 - AT" 1 = 0, 

H = vf+vl 

(39) 

Coincidentally, the equation for y is same as the equation for (D max /Dz) 2 (equation I21H . For a grid of values of [Di,T>2] I 
solve these equations for t]ez ,t)psf as a function of M and then numerically compute E(rj,r]EZ,TjpsF) and Ndet [A/\ X>i, O2]. 

I have carried out this numeric integration and then maximized Ndet with respect to N ■ I find that the optimum values 
of N vary from 0.54-0.60 for the entire range of ©1,2, with a mean of 0.57 and standard deviation of 0.02. So, in optimizing 
surveys I will just use the mean value of N = 0.57 since this only leads to a small error of at most 5%. The optimum value 
of Ndet varies by a factor of 2 (between the PSF and Exo-zodi limits, equation 1291 . so I have fit an equation to the numerical 
results: 



Ndet , 



1/2 , 1 1/2 



In fogf + 4a » 4(1 ~ ^ + ^sf(s 2 - 1)) - V% 



2(1~Vez+Vpsf(*o Z ~1)) 1 "" 

where ao = 0.4249 and a\ — 0.1287. This equation ranges from 12%/io (for PSF-dominated noise with t/psf = 1) to 24%/io 
(for Exo-zodi dominated noise with t\ez = 0) and it agrees with the numerically computed values of Ndet within < 1% 
fractional error! The number of stars being surveyed is N s = NN max , so the planet detection fraction of stars surveyed is 
actually [29%,44%,21%]/io in the [Zodi, Exo-zodi, PSF] cases; in practice the detection fraction will range between these values. 



9 OPTIMIZING THE CHOICE OF STARS 

Up until now, I have assumed that all stars are identical: same radius, same spectral shape, same luminosity, same planet 
frequency. However, surveys will be carried out for stars with a wide range of properties, so the question is: which are the 
best to observe? 

The answer to this should be obvious: prioritize the stars which have the highest probability for detection of planets. 
For a particular survey telescope, planet type, and exposure time, E(r], T>\, D2) should be computed for each star and then 
stars can be ranked in decreasing order of E. If one particular source of noise is expected to dominate for all stars, then 
computation of E can be sidestepped and stars can be ranked just based on the survey properties. 
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Figure 4. The dependence of the number of detected planets, Ndet on the number of stars surveyed, A/", both in terms of N max . Solid 
line is for the Zodi limit, dashed line is for the Exo-zodi limit, and dotted line is the PSF limit. 



To do this, I define a new parameter 

f3 XV - 1/s , (41) 

that depends solely on the properties of the stars at the observational wavelength used in computing r), 

ILq^D^t' 1 Zodi-limit 
P = { L^D~ 6 t ez Exo-zodi limit (42) 

L V D~ 6 PSF-limit 

Then, to choose stars I can simply rank stars from largest f3 to smallest /3 and observe in this order to maximize the number of 
detected planets. The parameter /3 has the advantage that it can be defined independent of the parameters of any particular 
survey so that the ranked list of stars simply depends on the wavelength of observation and dominant noise. 

For now, I will assume that the planet frequency is independent of stellar type and that stars are chosen of the same 
metallicity. For a given stellar type of uniform density one can show that N(> (3) = Npf3~ 3S . The normalization, Ng, is 
computed from 

{lI^Lq 3 /, 8 J drzp(Tz)r z 3 ^ s Zodi-limit 
Ll /2 J dT EZ p{T EZ )T- l z /2 Exo-zodi limit (43) 
Ll /2 PSF-limit, 

where p(tz)(1tz is the probability distributions of tz for the local zodiacal light; similarly, p(tez)<Itez is the probability 
distribution for exo-zodiacal light (both assumed to be independent of stellar spectral type) . I measure L v in units of 4-7T Jy 

pc 2 and D in units of pc. 

For computing N@ I use the present-day disc mass functio n of main s equen ce stars given bv iReid fc Hawlevl ll2005l) . 
the mass-temperature a nd radius-temperature relations given bv iHarmaned Jl98ctl . and stellar atmosphere models given by 
Irlauschildt et all 1^99!) • My computed Ng iX) is shown in Fig ure El where it is compared with a model of the disc luminosity 
func tion of dwarf stars in B, V, R, I and K Jjarrett et al.lll994) . I have checked that this result agrees with the sample of stars 
from lAllende Prieto &: Lambert! |[l999) (which is only complete for T e j j > 6000 K, so I have only compared in this range of 
temperatures). About half of Np is due to low temperature stars (T e // < 5000if) in the /-band, while more than one third 
is due to these stars in B and V. This indicates that late-type stars should be included in coronagraphic planet searches 
due to their greater number and thus greater proximity. This function can be used to compute the optimum wavelength of 
observation as shown below. 
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Figure 5. Normalization of j3 distribution as a function of wavelength for main sequence stars with f2 s = 47r, solid line is Zodi-limit 
while dashed-dot is Exo-zodi and PSF limits. The large dots are computed using the B, V, R, I and K band luminosity functions of 
nearby main sequence stars. For the Zodi and Exo-zodi limits, tz,ez has been set to one. The units of Na are chosen so that the specific 
luminosity has units of 47r Jy pc 2 and distance has units of pc. 



With N/3 in hand, the maximum number of stars to survey is 

8 / 3 n n 8 _t n ath 2 1 




3/11 



N° / - > GD? al T a [lAUY\- b @i£ A ) Zodi-limit 
N^gDtiT s \' 4 QJ^ A ) 1/[i Exo-zodi limit ( 44 ) 

N$gDl l T s \- 2 Co 1 \o 2 ®Tw A ) 1/3 PSF limit, 
where Q is defined in equation 1141 As before, if all three noise sources contribute then N max can be computed from equation 
1241 since T>\ , T>2 are weakly dependent on spectral type. Now, as in the case of a survey of stars uniform stellar type, the 
number of detected planets can be increased slightly if the number of stars surveyed is N = N a /N max =54-59%. So, the cutoff 
is determined by 

1 1/(3.5) 



3s 



(45) 



MNma 

So, all stars with /3 > /3 S can be chosen to optimize the survey. This guarantees that all stellar types will have a cutoff in 
distance at the same fraction of D max . Since I am assuming a constant exposure time per star, then the time will be divided 
up between various spectral types according to the integrand in equation 1431 For example, in the PSF-limited case, stars of 

1/2 

luminosity L„ will have time spent dT s oc (dn/dM)Lj dM. 



10 OPTIMIZING THE WAVELENGTH OF OBSERVATION 

Since A^et is a function of wavelength, we need to compute the wavelength dependence of Ndst/N rnax and N max . To compute 
the first quantity, T>\ and T>2 should be computed from equation 1231 for G-type stars (which are most common) at a range 
of wavelengths. Since Z?i,2 are weakly dependent on spectral type, these can be used for optimizing a survey of all stellar 
types. Then, t]ez and t)psf can be computed from equations 1391 by computing R from the fifth equation, solving for y from 
the fourth equation (a quartic), and then plugging y into the first two equations. Then, equations 1251 and 1441 can be used 
to compute N max from the other telescope, stellar and survey parameters, while equation 1401 can be used to compute the 
number of detected planets from t)ez,psf and N max . This procedure yields Ndet{\), which I carry out for a range of satellites 
discussed in the next section. 
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The general wavelength scaling can be derived if one of the three noise limits dominates since N max with wavelength as 

{(ivj /3 eA- 6 p!) 3/11 Zodi limit 
(N^eX-^l) 1 ^ Exo-zodi limit (46) 
(Nje\- 2 p 2 x ) llS PSF limit 

assuming A In A = constant. I plot N max versus wavelen gth for Earth a nd Jupiter albedos shown in Figure |S| The albedo of 
the Earth I take from obser vations of earthshine bv lWoolf et alJ (120021) and at shorte r wavelengths from direct observations 
of the Earth with G OME teurrows et aljll999h . The albedo of Jupiter I take from iKarkoschkal (Il994h in the optical and 
IWaeener et alJ l)l985h in the ultraviolet. I have assumed a 20% bandpass for computing these curves and I assume that e\ is 
independent of wavelength. 

Figure |2| shows that the peak of Np is near 1.6 micron (this is due to the large number of stars cooler than the Sun). 
However, due to the A - ' 6,4,2 ' term, the peak in the number of detected planets is at a much shorter wavelength, A./ = 
[3961, 4785, 5826]A for Jovian albedo planets and A® = [3953, 4615, 4794]A for Earth albedo planets in the [Zodi,Ezo-zodi,PSF] 
noise limits, respectively. For Earth-like planets a sharp drop occurs shortward of 3000 A due to the absorption opacity of 
ozone. Thus, surveys in the B and V bands should have the highest efficiency in detecting planets. 



11 COMPARING PROPOSED SATELLITES 

With the analytic relations derived in the previous sections in hand, I can now compare various proposed coronagraphic planet 
imaging satellites. A survey of the SPIE conference proceedings "Future EUV/UV and Visible Space Astrophysics Missions 
and Instrumentation" and "High-Contrast Im aging for E xo-Planet Detection" turns up at lea st six pla nned coronagraphic 
satellites: EPIC/OP D jMennesson et alJ | 2003h ECLIPSE jTraueer et al.ll2003l: IHuII et aljEool. TPF-C I Brown et al]l2003t 
iTraub fc etafl boOfft ESPI l|Lvon et alJ booA . UMBRAs Jschultz fc et alJbooA ExPO llGezari et al.ll2003h . and TOPS 
l|Giivmi^^t^ni2006rl . For comparison, I also include the Hubble Space Telescope lBrow n~fc^um?wa^990F) and the Hubble 
Space Telescope with a corrective secondary mirror, HST* as proposed bv lMalbet et alji^95TT "l also include the "DREAM" 
satellite which is my version of what an ideal design for TPF-C would be capable of. The parameters of these different missions 
reported in these citations are listed in Tabled 

To compare these various satellites, I def ine a figure of merit which I call "planet imaging power" , or PIP, (analogous to 
the "A-fl" etendue defined for large surveys, Iciaver et alJl20o3 : lKaiseill2004F ) . It is given simply by 

f(A 2 eV 3 e7vl?A) 1/11 Zodi limit, 

PIP = I (Dl%e A Q^ A fl 9 Exo-Zodi limit, (47) 

[(Dl^e^C^ 9 PSF limit. 
The number of planets that may be detected with a given telescope is simply proportional to PIP assuming that a particular 
source of noise dominates. I have computed these for each proposed satellite in Table Q which shows that TPF-C is the 
top-ranked proposed satellite in all three noise limits. 

To compare these satellites, I assume a throughput of e —3% (independent of wavelength - although for TOPS and 
DREAM I assume an efficiency of 20%), a mission duration of T a = 3yr, a solid angle of f2 s = 4-7T, and a planet albedo 
identical to either the Earth or Jupiter. I compare four surveys (1) a 3-year survey of solely G-type (Sun-like) stars optimized 
for either jupiters or earths; (2) a 3-year survey of all nearby main-sequence stars optimized for either jupiters or earths. 
Additional assumed parameters are A In A = 0.2, (S/N)det = 10, tz = 6.8 x 10~ 10 sr" 1 (which is the weighted value in the 
anti-solar direction), tez = 2.7 x 10~ 9 sr" 1 (which assumes an edge-on zodiacal cloud like the Sun's around every star), 
a quasi-Lambert phase function, and, for the G-star survey, I assume a stellar number density of n* = 5 x 10 -3 pc -3 . For 
Jupiter-albedo planets I assume R p — 7.15 x 10 9 cm and and for Earth-albedo planets I assume R p — 6.38 x 10 s cm. I compute 
the optimum T exp , N max , Ndet, and /<j ei = Ndet /(N s f io) for each proposed survey, listed in Tables I2T51 (max is abbreviated 
m for some quantities in the table) . Two examples for the G-dwarf survey of earth-sized planets with TPF and a survey of 
all main sequence stars for jupiter-sized planets with TOPS are shown in Figure El For the computed exposure times, I also 
compute D max and r max for Sun-like stars. I assume that there is one planet per decade of radius for every star (that is, 
/io = 1). Since the current data favors /io ~ 10% for Jovian planets (Tables [5] and these numbers should be reduced by a 
factor of ~ /io. 

Since each telescope has slightly different design parameters, the values in Tables 12151 are only an approximate guide. 
However, a few points are striking. All of the detection efficiencies are between 21-32% of the average number of stars per 
decade of radius for the surveyed stars. Although the ECLIPSE and EPIC telescopes have similar sizes, the EPIC telescope 
is much more powerful due to its smaller inner working angle. For finding terrestrial planets about G dwarfs, the TPF-C is 
the only one with a statistically significant chance due to its large aperture size. More recent designs for TPF-C include a 
larger mirror than I have used here, so this number may be an underestimate. With a 4m mirror, an optimized survey of G 
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Figure 6. The maximum number of survey stars, N max , and number of detected planets, N aet , as a function of wavelength for a 
population of planets with identical albedos to Jupiter (top) and Earth (bottom). The solid line is N max for the Zodi-limited noise case, 
the dashed line is the Exo-zodi limit, and the dotted line is the PSF-limit. The dash-dot line is Nmax from the combined noise limits 
for a survey with parameters of the TOPS telescope surveying all main-sequence stars (top panel) and parameters of TPF- C surveying 
only G dwarfs (bottom panel). The dark solid lines are the number of detected planets, N^ e t, for an optimized survey with /io = 1. 



dwarfs with TPF-C is most sensitive to planets at ~2 AU, outside the habitable zone (if fewer stars are surveyed, then TPF-C 
can be optimized for the habitable zone, but many fewer planets will be detected). Jupiter-sized planets are optimized for 
detection at a few AU, but will be detected at a large range of radii. Comparing tables 0] and shows that there is a factor of 
~3-5 increase in the number of detected earth-sized planets if all stars surveyed rather than just G-dwarfs. It is remarkable 
that the much smaller 1.2m TOPS telescope will have a significant chance of detecting a few Earth-sized planets if /io is 
not too small. The optimum wavelengths of detection range from 4000-4900 A for earth-albedo planets and 4400-5900 A for 
Jupiter-albedo planets; thus, the B and V bands are optimal for planet detection. The last two columns give the ratios of the 
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Table 1. Proposed satellite parameters. 



Satellite Q IWA D tel C PIP Z PIPez PIPpsf 
(m) @550nm 



ECLIPSE 


4.0 


1.8 


l.e- 


-09 


0.3 


0.4 


0.16 


EPIC 


1.5 


1.5 


l.c- 


-09 


0.6 


0.5 


0.41 


ESPI 


4.0 


1.5 


3.e- 


-07 


0.2 


0.3 


0.02 


ExPO 


2.9 


3.0 


l.e-09 


1.3 


1.1 


0.68 


TPF-C 


3.9 


4.0 


5.e- 


-11 


1.8 


2.0 


2.20 


UMBRAS 


3.5 


1.0 


l.c- 


-OS 


0.1 


0.1 


0.03 


HST 


3.0 


2.4 


l.c- 


-06 


0.8 


0.9 


0.04 


HST* 


3.0 


2.4 


5.e- 


-10 


0.8 


0.9 


0.52 


TOPS 


1.5 


1.2 


l.e- 


■10 


0.6 


0.6 


1.06 


DREAM 


1.5 


S.O 


5.e- 


-11 


38.7 


28.6 


59.17 



Table 2. Expected number of detectable Jovian planets surveying all stars. 



Satellite 




N a 


Ndetl 7l0 


fdet 


Am 


Dm 






T>2 


Acronym 


(hr) 






(%) 


(A) 


(pc) 


(AU) 


















m 


no 






Equation(s) 


I20I24I44I 






mum® 


ESIiSlifll 


ECLIPSE 


41.7 


630 


152 


24.2 


5189.3 


23.6 


6.5 


0.69 


0.94 


EPIC 


18.9 


1392 


366 


26.4 


4897.5 


29.7 


3.5 


0.85 


0.86 


ESPI 


311.8 


84 


17 


20.9 


5848.7 


12.3 


4.6 


0.62 


2.11 


ExPO 


9.7 


2710 


662 


24.4 


5189.3 


38.0 


4.5 


0.77 


0.95 


TPF-C 


5.5 


4790 


1368 


28.6 


4412.9 


44.6 


1.6 


0.79 


0.68 


UMBRAS 


198.9 


132 


28 


21.7 


5584.0 


14.2 


6.6 


0.65 


1.22 


HST 


123.9 


212 


44 


20.9 


5848.7 


16.7 


2.9 


0.68 


2.56 


HST* 


14.5 


1810 


465 


25.7 


4897.5 


33.0 


4.8 


0.76 


0.87 


TOPS 


15.5 


1691 


503 


29.8 


4278.7 


31.0 


3.9 


0.82 


0.59 


DREAM 


0.3 


105044 


32520 


31.0 


4412.9 


117.9 


2.3 


0.97 


0.61 



distance maxima in the different noise limits. For all of the surveys the different sources of noise contribute at a similar level - 
this is not surprising as the telescopes were designed with the realization that it does not pay to reduce the PSF wings much 
below the Zodi and Exo-zodi limits. This is apparent in the both panels of Figure ©where the Zodi and PSF noise contribute 
equally for the TOPS survey of jupiters, while Zodi and Exo-zodi noise contribute equally for the TPF-C survey of earths 
around G-dwarfs. 



12 OPTIMIZING THE METALLICITY OF SURVEY STARS 

Since it is known that planets of period less than 4 ye ars show a frequency proportional to the square of the metallicity of the 
host star ^Fischer fc Valentill200,4 ISantos et al.ll2004l ). the question arises as to how to distribute observing time among stars 



Table 3. Expected number of detectable Jovian planets around Sun-like stars. 



Satellite 


Texp 


N a 


AW/io 


fdet 


A ,77 


Dm 




v 1 


V 2 


Acronym 


(hr) 






(%) 


(A) 


(PC) 


(AU) 








HI 




1401391 






m 


1271 






Equation(s) 


IMI1GI44I 








ESIifllifll 


ECLIPSE 


130.8 


200 


49 


24.8 


5149.4 


27.7 


7.6 


0.66 


0.90 


EPIC 


63.5 


413 


110 


26.8 


4897.5 


35.3 


4.1 


0.81 


0.81 


ESPI 


911.4 


28 


6 


20.9 


5562.5 


14.7 


5.2 


0.60 


2.07 


ExPO 


30.9 


850 


211 


24.9 


5149.4 


44.9 


5.3 


0.74 


0.91 


TPF-C 


19.0 


1382 


399 


28.9 


4568.8 


52.6 


5.6 


0.74 


0.64 


UMBRAS 


592.4 


44 


9 


21.9 


5562.5 


16.9 


7.8 


0.62 


1.17 


HST 


362.1 


72 


15 


20.9 


5562.5 


20.1 


3.3 


0.65 


2.51 


HST* 


47.3 


555 


115 


26.2 


4897.5 


38.9 


5.7 


0.73 


0.82 


TOPS 


55.2 


476 


141 


29.8 


4533.7 


36.9 


5.0 


0.77 


0.54 


DREAM 


1.0 


27282 


8365 


30.7 


4568.8 


142.8 


2.9 


0.90 


0.57 
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Table 4. Expected number of detectable Earth-like planets surveying all stars. 



Satellite 
Acronym 


(day) 


N a 


Ndet / flO 


fdet 
(%) 


Am 

(A) 


D m 

(PC) 


(AU) 


Pi 


V 2 




m 




I4QI39I 






TM 


1551 






Equation(s) 


1201^41441 






1231141461 


I23li3li8l 


Hi O Hi 1 JT O Hi 


32 1 




7 


22 


AQfi'i O 


q n 


9 3 
Z .o 


n qi 


1 97 


Hirl O 


la. 1 


o.> 




OA 

Z4.o 


4ool.o 


11. o 


l.o 


ill 
i.ii 


i n 
1 . iz 


ESPI 


299.6 


3 





20.9 


4903.2 


4.4 


1.4 


0.83 


2.99 


ExPO 


7.5 


145 


32 


22.5 


4903.2 


14.4 


1.6 


1.02 


1.28 


TPF-C 


3.0 


370 


101 


27.4 


4412.6 


18.2 


1.9 


1.02 


0.89 


UMBRAS 


182.2 


6 


1 


21.0 


4903.2 


5.2 


2.1 


0.86 


1.69 


HST 


119.1 


9 


1 


20.9 


4903.2 


6.0 


0.9 


0.90 


3.62 


HST* 


10.2 


107 


25 


23.6 


4881.3 


12.8 


1.9 


1.00 


1.13 


TOPS 


7.5 


145 


43 


30.0 


4079.4 


12.5 


1.5 


1.09 


0.79 


DREAM 


0.1 


8634 


2768 


32.1 


4181.3 


AAA 


0.8 


1.28 


0.83 



Table 5. Expected number of detectable Earth-like planets surveying Sun-like stars. 



Satellite 




N a 


Ndet/ f 10 


fdet 


Am 


D m 




Pi 


v 2 


Acronym 


(day) 






(%) 


(A) 


(pc) 


(AU) 








El 




1401391 






VM 


na 






Equation(s) 


12(11191441 






ESIiSlifll 


I2SI1SI1SI 


ECLIPSE 


94.4 


11 


2 


21.7 


4522.8 


10.0 


2.4 


0.96 


1.38 


EPIC 


43.6 


25 


6 


24.2 


4422.5 


13.0 


1.4 


1.18 


1.23 


ESPI 


848.6 


1 





20.9 


4533.0 


4.8 


1.4 


0.88 


3.26 


ExPO 


22.7 


18 


10 


22.2 


4522.8 


16.1 


1.7 


1.08 


1.39 


TPF-C 


11.0 


99 


27 


27.3 


4422.5 


20.8 


2.1 


1.06 


0.93 


UMBRAS 


518.8 


2 





21.0 


4533.0 


5.6 


2.1 


0.92 


1.85 


HST 


337.4 


3 





20.9 


4533.0 


6.5 


0.9 


0.96 


3.94 


HST* 


31.7 


34 


7 


23.0 


4422.5 


14.5 


1.9 


1.06 


1.24 


TOPS 


31.8 


34 


10 


30.3 


4422.5 


14.6 


1.9 


1.10 


0.78 


DREAM 


0.6 


1718 


556 


32.4 


4422.5 


52.9 


1.0 


1.30 


0.82 



of different metallicities. Although coronagraphic surveys for giant planets may tend to find planets at periods longer than a 
few years, it may be reasonable to assume that the planet frequency-metallicity correlation will also hold for the distribution 
of planets discovered in a coronagraphic survey. I will assume that the exposure time is adjusted in proportion to the expected 
number of planets per star raised to some power 7, T exp (X) oc /7 (X =[Fe/H]). If I assume that the magnitude of a given 
type of star is independent of the metallicity, then the number of planets found about that stellar type is 

r rl+3-yS/(l+3S) dn, ,y 

Ndet OC J%+™ ■ (48) 

U ho-dx dJi ) 

By differentiating with respect 7, the number of planets is maximized for 7=1 + 35, mea ning that each X bin sho uld be 
weighted by metallicity to the ~third power. Using the metallicity distribution measured bv lValenti fc Fischerl I^OOfil) I find 
that the number of planets can be increased by 14% or 19% for 7 = 1 + 35 (with 8 = 1/8 or 1/6) versus 7 = (i.e. no 
preference for higher metallicity stars). 



13 DISCUSSION AND CONCLUSIONS 

Since detection of planets is the first-order goal of a coronagraphic survey, maximizing Ndet has some interesting implications: 
1) Survey duration: The total number of planets detected scales as Ndet oc T^ 11 in the Zodi limit or oc T s 1/3 in the 
Exo-zodi and PSF limits, which means that about ten times as much time must be spent to double the number of detected 
planets. 

3) Telescope design optimization: The number of planets detected scales as Ndet oc PIP (given in eauation l47|l in the three 
different noise limits. The greatest gains come from first increasing the telescope aperture, second decreasing the inner working 
angle, and third decreasing the PSF wing contrast (if the PSF limits the detection) or increasing the telescope efficiency. 

4) Optimum wavelength: This must be computed by combining the optimum number of stars to survey in all three noise 
limits. For Earth albedo the peak tends to lie in the B band, while for Jupiter albedo it lies in the V band (however, different 
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telescope and survey parameters might shift the optimum wavelength to another region). In all cases the optimum wavelength 
is shortward of the peak of the stellar photon number flux. 

5) Stellar spectral types: The number of planets detected scales with spectral type as Ndet oc dnV 11 L^ 11 in the Zodi noise 
dominated limit, while Ndet oc dn'J i L\!' i in the Exo-zodi and PSF dominated limits. Stars of G type dominate the number 
of detected planets (assuming that the frequency of planets is independent of stellar type). If main sequence stars of all types 
are surveyed, then the number of planets can be increased by a factor of ~ 3 over surveys that only target G-type stars. 

6) Planet size: I have computed the number of detected planets assuming a single planet size (equation However, we 
know that there will be a distribution of sizes of planets, in which case the number of planets detected at other sizes scales 
as ^| - -Rp 3/ ' 2 2 ' f° r smaller planets (in the [Zodi,Exo-zodi/PSF] limits), while for larger planets it has a more complicated 
scaling. So, for example, Earth-sized planets would need to be ~10 2 times more abundant than giants to be detected in 
similar quantities by a survey optimized for giants. 

7) Signal-to-noise: Since the number of planets in an optimized survey above a signal-to-noise (S/N)o scales as the volume 
within which those planets can be detected, then 



N det [S/N > (S/N) ] = N det [S/N > (S/N) det ] { V CfWo J ' (49) 



3/4 

Zodi-limit, 



(S/N) d 
(S/N) 



Exo-zodi/PSF limits. 



8) Planet-star distance: Since D max sets the maximum distance to which planets may be detected, r max = 6iwAD ma x/ sina,, 
sets the planet-star separation at which the bulk of the planets will be detected (most planets are detected close to r ma x since 
this is where the volume of the survey is maximized). Since G type stars dominate the number of detections, r max should 
be computed for G stars in all three noise limits, and the minimum value should provide a good estimate of the peak of the 
planet-star distance distribution. If there is a reason to believe that the target planet population is not distributed uniformly in 
log(r), then the optimization of the observing strategy will change only slightly. I have run some test cases with a planet-star 
separation distribution that is a power law with radius, and I find that the optimum N is still around 0.6. The reason for 
this is that the volume term dominates - at large distances the number of stars surveyed grows rapidly, so it does not pay 
to greatly decrease the number of surveyed stars. However, the values of N det can change significantly for a non log-uniform 
distribution of planet-star separations. For a planet-star separation that scales as df/dr oc r~^ 1+ ^ , total number of planets 
detected will scale as r^ ax , and so the specific value of r max for a given telescope and survey strategy will determine how 
many planets are detected. 

If one is concerned with targeting habitable-zone planets, then a different optimization s trategy is req uired (which will 
likely involve observing more nearby stars); I will defer a study of this to future work (see iBrownl [20051 . for a numerical 
approach). 

9) Fraction of detectable planets: I computed that if all stars are surveyed within D max , then ~ 21 — 44% of the number 
of planets per decade in radius will be detected (even though a smaller number of stars are surveyed than in an unoptimized 
survey ) ; the satellites modeled here show a range of 21-32%. Multiple visits might increase this yield as discussed bv IBrownl 

(2004) , which likely requires a Monte-Carlo approach, so I defer a study of this issue to future work. 

10) Metallicity: If the planet frequency increases as the square of metallicity (as it does for the known extrasolar giant 
planets), then the exposure time should be scaled as the ~cube of the metallicity. This increases the number of detected 
planets by about 15-20% over the number detected with a uniform exposure time. 

11) Proposed telescopes: The TPF-C is by far the most powerful and ambitious of the proposed coronagraphic imaging 
telescopes. However, a much smaller telescope TOPS has a significant chance of detecting Earth-size/albedo planets if they 
are not too rare. The main competitive strengths of TOPS axe (a) smaller inner working angle and (b) high throughput. Thus, 
a modified version of TPF-C with a larger mirror, higher throughput, and smaller inner working angle (I dub "DREAM") 
would allow more than an order of magnitude increase in the number of detected Earth-like planets and would have a peak 
sensitivity in the habitable zone. 

Finally, I wish to comment on the few assumptions I have made which should be studied in future work. I have assumed that 
the albedo and phase function is independent of the separation between planet and star and spectral type of t he host star, and 
that t he phase function is independent of wavelength. The simulated Jupiter-mass planetary atmospheres in lSudarskv et, all 

( 2005) show an inverse square-law dependence with radius in the V-band from ~ 2 — 20 AU, show a nearly constant phase 
function at > 45° from 0.2-4 AU, and show a geometric albedo that differs by at most 0.2 in the optical over 2-15 AU from 
the albedo at 6 AU. In regions that differ from an inverse square law, the results in this paper will change quantitatively, 
but probably not qualitatively. The Lambert phase function shape is a ve ry poor description of cloudy atmospheres which 
can have a peak in the brightness at crescent phase l)Sudarskv et alJl2005h : however, this effect is more significant at longer 
wavelengths which are unfavorable due to poorer angular resolution. 

I have assumed that the PSF has a hard "edge" at the inner working angle within which no planet can be detected, 
and that the wings of the PSF have a constant contrast. This is also clearly an oversimplification, but a full study will 
require knowing precisely the properties of a given instrument. I have assumed that the telescope is circular, while plans for 
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TPF-C indicate that an oblong telescope may give better detection properties without much added cost. My analysis can 
accommodate an oblong mirror which has an inner working angle set by the long axis, which can be taken to be 2R te i- Then, 
the area of the telescope is smaller than nR 2 e i by the ratio of the axis lengths which can be included in the efficiency factor, e. 
Observations at several roll angles can be used to circularize the inner working angle of an oblong telescope which also leads 
to a hit in the efficiency by the inverse of the number of roll angles. 

I have assumed that the telescope efficiency is independent o f wavelength. This m ight be achieved, for example, by 
observations with a superconducting tunneling junction (ST J) array JPeacock et alJl998f) . but until then one must repeat my 
analysis with the efficiency of a particular telescope folded in. STJ detectors have the additional advantage that they measure 
the energy of each photon allowing detection of broad spectral features and subtraction of speckles which have different widths 
in each waveband. I have also assumed that there is no limit on exposure time - it may be that pointing control, cosmic rays, 
or other factors will limit the exposure time per star, which has not been taken into account in this analysis. 

I have assumed that the planet frequency is independent of spectral type; however, enough data to quant ify the variation 
of planet frequency with spectral type currently do not exist, although some predictions have been made ijlda fc Linll2004l : 
lLaughlin et al-lfeoO^I) . One of the most significant uncertainties in these calculations is the strength of the Exo-zodiacal light, 
of which there are currently no observational constraints, ft is likely that the Exo-zodiacal dust distribution will depend on 
the age, spectral type, planet and minor planet distributions around given stars, so there is no simple way to predict or 
parameterize this source of noise. Further studies in the mid-infrared would be very helpful in constraining this factor and 
assessing how important it will be in limiting planet detection. 

f have ignored the effect of Oowa since it only causes obscuration of planets orbiting the nearest stars, blocking a fraction 
~ {Oiw a/8ow a) 2 which is typically quite a small number. 1 have primarily concentrated on a constant exposure time for all 
stars, but weighting the exposure time towards nearer stars can cause a slight increase in the signal-to-noise for each detection 
and increase the number of nearby planets detected at close distances to the host star around stars closer to the Sun without 
decreasing the number of detected planets significantly. 

I have assumed for the blind surveys that there is no prior information about the existence of planets about these stars. 
However, given that of order one thousand stars have been surveyed to date with the radial velocity technique, it is likely 
that some parameter space for planets can be excluded for som e survey stars. Currently the surveys are complete to about 
K ~ f 5 m/s for long period planets, ~ 3 years while the peak of the surveys for Giant planets is about 3-f 5 

years, so the current constraints will likely rule out only a small region of parameter space. However, as radial velocity surveys 
improve their sensitivity and time baseline these constraints will become more stringent. Incorporating this information into 
target selection and exposure time will be left for future work. 

If you wish to use the codes used to compute the figures or tables, please visit the author's web page at the University 
of Washington Astronomy Department. 



14 APPENDIX 

To make the notation in this paper easier to navigate, here are two tables 16171 referencing the symbols used throughout, the 
units of each quantity, and the section or equation the first time each symbol appears in the paper. 
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Table 6. Definition of Roman symbols. 



oymbol 


U SagC 


Units 


Ej C|U. cit 1 H 


a-0, a l 


Numbers used in approximation to N det /N max 


- 


13771 


C(X) 


Intensity contrast of PSF 


— 





Co 


Intensity contrast of PSF at a fiducial wavelength 


— 


i 


D 


Distance of star from the Earth 


cm or pc 


in 


£>i 


Ratio of Dz to Dez 


— 


EH 




Ratio of Dz to Dpsf 


- 


l23l 


Dmax 


Maximum distance a planet can be detected 


cm or pc 


U 


D s 


Maximum distance of survey 


cm or pc 


El 


O z 


Dmax in the Zodi noise dominated limit 


cm or pc 


IT31 


D EZ 


Dmax in the Exo-zodi noise dominated limit 


cm or pc 


rm 


DpSF 


Dm a x in the PSF noise dominated limit 


cm or pc 


rm 


D 


Distance of star from the Earth 


cm or pc 


□ 


E(V) 


Expection value of the number of planets that can be detected 


- 


vm 


E 


Approximation to E(rj) 


- 


rm 


ho 


Expected number of planet per decade of r 


_ 


□ 


fdet 


Fraction of stars observed with detected planets, per decade radius 


— 


JED 


G 


Auxiliary quantity used in defining D ma x 


_ 


I22I 


Q 


Quantity used in computing D max and N ma x 


cm 4 erg -1 s -1 


rm 


h 


Planck's constant 


erg s 


□ 


H 


Heaviside step function 


_ 


1 


L u 


Star's specific luminosity 


erg s — 1 Hz -1 


- 


L v ,® 


Sun's specific luminosity 


erg s — 1 Hz -1 




M 


Stellar mass 


M 


□ 


jV 


Ratio of N a to N max 




rem 


Ndet 


Number of detected planets 


_ 


in 


Nmax 


Maximum number of stars that can be surveyed 


_ 




Nz,EZ,PSF 


Maximum number of stars that can be surveyed in 3 noise limits 




1191441 


N s 


Number of stars in survey 




□ 


N 


Normalization of f) distribution 




PHil 


n* 


Number density of stars of a particular spectral type 


pc -3 


JO 


Pa 


Geometric albedo of planet as a function of wavelength 


— 


m 


PIP 


Planet imaging power (in each noise limit) 


various 


E71 


Qez 


Variance of noise from exo-zodiacal light 


_ 


m 


Qz 


Variance of noise from Zodiacal light 


- 


EH 


Qpsf 


Variance of noise from stellar PSF 


_ 


ED 


Qb 


Variance of noise from other backgrounds 


_ 


6 


Qp 


Number of photons detected from planet in an exposure 


— 


m 


Q* 


Number of photons detected from star in an exposure 


- 


m 


r 


Planet-star separation 


cm or AU 


m 


n 


Minimum r at which 8 = Bjwa 


cm or AU 


rm 


Tl 


Maximum r at which S/N = S/N det 


cm or AU 


1521 


T2,Z 


Maximum r at which S/N = S/N det in Zodi noise limit 


cm or AU 


rm 


T2,EZ 


Maximum r at which S/N = S/N det in Exo-zodi noise limit 


cm or AU 


rm 


r2,PSF 


Maximum r at which S/N = S/Ndet in PSF noise limit 


cm or AU 


rm 


Tmax 


Planet-star separation at maximum elongation 


cm or AU 


rm 


K 


Auxiliary quantity used in defining r\EZi VPSF 


- 


l39l 


Hp 


Planet radius 


cm 


□ 


Rtel 


Radius of telescope 


cm 


□ 


S 


Quantity used in defining E(rf) 


_ 


1251 


S/N 


Signal-to-noise ratio 


- 


□ 


S/N 


Fiducial signal-to-noise ratio 


- 


jm 


S/N det 


S/N required for detection 




m 


Sfac 


PSF sharpness 




8 


Texp 


Exposure time per star 


sec 


H 


T s 


Total time of survey 


sec 


JU 


V 


Auxiliary quantity used in defining D ma x 




I22I 


w 


Auxiliary quantity used in defining dN^/dQ in PSF noise limit 




rm 


x± 


Quantity used in defining E{rf) 




l28l 


X 


[Fe/H] log of metallicity 


dex 


□ 
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Table 7. Definition of Greek symbols. 



O V lilUUl 


T T G fro 


Units 


TV,/~i iiof i / \ ri 

±_J*J LlcilrlWll 


a 


Planet phase angle 


rad 


□ 


a± 


Phase angle where resolution and flux limits cross 


rad 


1551 




Phase angle where planet-star sky angle is maximum 


rad 




P 


Parameter used to rank priority of observing different stars in each noise limit 


various 


l42l 


Ps 


Value of P at D s 


_ 


1351 


7 


Exponent of exposure time versus planet frequency as a function of metallicity 


_ 


1121 


r 


Quantity used in defining E(r]) 




1551 


Six) 


Dirac delta function 


— 


□ 


8 


Exponent of T exp in equations for Dmax 


- 


1T51 


A 


Quantity used in defining E(rf) 


- 


1551 


A In A 


Fractional bandwidth of telescope 


- 


El 


e 


Total throughput 







c 


Ratio of r to r max 






'/ 


Ratio of D to D m ax 


_ 


El 


Vz 


Ratio of D m ax to D z 


- 


1501 


VEZ 


Ratio of D m ax to D E z 


- 


1501 


riPSF 


Ratio of Dmax to D PSF 


- 


1501 


e 


Sky angle separation of planet & star 


rad 


El 




Inner working angle of coronagraph 


rad 


1 


®IWA 


Dimensionless inner working angle 


- 


m 


K 


Power law exponent for exposure time as a function of distance 




B 


A 


Wavelength 


A 


3 


Ao 


Reference wavelength for C(A) 




EI 


S 


Auxiliary quantity used in defining dNdet/dC m PSF noise limit 




1551 




Power law exponent of df /d In r versus r 




EI 


TZ 


Relates solar flux to zodiacal surface brightness 


sr _1 


m 


T~EZ 


Relates stellar flux at tangent to line of sight to exo-zodiacal surface brightness 


sr" 1 


m 




Quantity used in defining E(rf) 




1551 


*(a) 


Planet phase function 




m 




Solid angle of planet survey 


rad 


□ 
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